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If dark energy interacts with dark matter, this gives a new approach to the coincidence problem. 
But interacting dark energy models can suffer from pathologies. We consider the case where the 
dark energy is modelled as a fluid with constant equation of state parameter w. Non-interacting 
constant-w models are well behaved in the background and in the perturbed universe. But the 
combination of constant w and a simple interaction with dark matter leads to an instability in the 
dark sector perturbations at early times: the curvature perturbation blows up on super-Hubble 
scales. Our results underline how important it is to carefully analyze the relativistic perturbations 
when considering models of coupled dark energy. The instability that we find has been missed in 
some previous work where the perturbations were not consistently treated. The unstable mode 
dominates even if adiabatic initial conditions are used. The instability also arises regardless of how 
weak the coupling is. This non-adiabatic instability is different from previously discovered adiabatic 
^—2 . instabilities on small scales in the strong-coupling regime. 

^ ' PACS numbers: 95.36,+x, 98.70.Vc, 98.80.Cq 
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O ■ I. INTRODUCTION 

<^ ! 

i— i, In the standard cosmological model, dark energy and dark matter are the dominant sources in the evolution of 

i-P ■ the late universe. They are currently only indirectly detected via their gravitational effects, and this produces an 

y^l important degeneracy [l|. In particular, there could be a coupling between dark energy and dark matter without 

Q i violating observational constraints. A coupling in the dark sector could help to explain why the dark energy only 

comes to dominate after galaxy formation. But some of these models may be ruled out by instabilities that are not 

CZ2 , apparent in the background solution. 

i i ' Various forms of coupling have been considered (see e.g. [3, y, 0| and references therein). A general coupling may 

be described in the background by the energy balance equations of cold dark matter (c) and dark energy (x) , 

£>. . p'c = -3Hp c - aQ , (1) 

(N ! p' x = -m(l + w x )p x + aQ, (2) 

fvq , where w x — P x j p x , H. = dlna/dr and r is conformal time, with ds 2 = a 2 (—dr 2 + dx 2 ). Here Q is the rate of energy 
density transfer, so that Q > (< 0) implies that the direction of energy transfer is dark matter — > dark energy (dark 
i* \ energy — > dark matter). 

The density evolution in the dark sector deviates from the standard case. We can use effective equation of state 
parameters for the dark sector to describe the equivalent uncoupled model in the background: writing p' c + 37Y(1 + 
O ■ Wc,cft)Pc = and p' x + 37Y(1 + w x , cS )p x = 0, we have 

aQ aQ 

W c ,cS = 77T- , Wx,e& = W x - — . (3) 

otipc onp x 

When Q > 0, we have w CtC g > 0, so that dark matter redshifts faster than a" 3 , while W x , e s < w x , so that dark energy 
has more accelerating power. The opposite holds for Q < 0. When Q > 0, the coupled dark energy can behave like an 
uncoupled "phantom" model, i.e., with w x ^r < — 1, but without the usual problems associated with phantom dark 
energy 0]. 

In order to avoid stringent "fifth-force" constraints, we assume that baryons (b) and photons (7) are not coupled 
to dark energy and are separately conserved, and we assume the same for neutrinos (y). So the balance equation for 
fluid A is 

p' A = -3W(1 + w A )p A + clQa , (4) 

with Qt, = Q 7 = Q u = and Q c — —Q — —Q x =/= 0. The Friedman equation is 

o o7TGr 2 / \ / r \ 

n = ^— a (P~t + Pk + Pb + Pc + Px)- (5) 

Once a form for Q is given, the background dynamics are fully determined by the above equations, and typically 
the analysis focuses on the possibility of accelerating attractor solutions (for recent work with further references, see 
e -g- @iS|)- The models may also be tested against geometric observational constraints (see e.g. [2J]). 

In the perturbed universe, there are subtleties and complications that do not arise for the background dynamics. 



X 



• Firstly, one needs a covariant form for the dark sector energy-momentum transfer that holds in an inhomogeneous 
universe, and reduces to the background form in a Friedman-Robertson- Walker (FRW) universe. For example, 
if one uses the ansatz Q = Qoa n , then the background dynamics can be determined and the parameters Qo,n 
can be constrained by geometric observations. However, there is no covariant form for such ad hoc ansatzes, and 
therefore one is unable to compute the perturbations - no consistent cosmological model can be constructed on 
the basis of such ansatzes. 

• Secondly, one needs to ensure that dark energy perturbations are stable, i.e., <r sx > where c sx is the dark 
energy sound speed (the speed at which fluctuations propagate) . For a scalar field model of dark energy, c 2 x = 1 
follows without assumptions 6] . But for fluid models as used here, we need to impose c 2 x > by hand, so that 
the dark energy fluid is effectively non-adiabatic. 

The sound speed problem applies equally to uncoupled dark energy, but since the coupling itself can introduce 
non-adiabatic modes, the issue is even more important in the coupled case. 

Here we consider a dark energy fluid, with w x = const, that is coupled to dark matter via a covariant energy- 
momentum transfer four- vector Q^, which reduces in the background to ±o,~ 1 QSq, where Q is a simple function of 
energy density. We show that the gauge-invariant curvature perturbation has a super- Hubble instability in the early 
radiation era, no matter how small the coupling is, and even if adiabatic initial conditions are used. This rules out 
these models. It appears that constant-u^ fluid models of dark energy, even with the imposition of c 2 x = 1, are 
unstable to couplings with the dark matter. The non-adiabatic large-scale instability that we find is different from 
the small-scale instabilities in the strong coupling and adiabatic regime that have been previously discussed pj Q • 

In order to avoid the large-scale instability, w x must increase sufficiently in the early radiation and matter eras. In 
other words, the simple constant-u^ fluid model cannot be extended to the early radiation era. This is in contrast 
to the case of uncoupled dark energy, where constant-it^ fluid models are well behaved in the primordial universe. 
A dynamically evolving quintessence field will likely avoid the instability we find, since w x typically does not remain 
constant back to the early radiation era. 

The plan of the paper is as follows. In Sec. II, we present the density and velocity perturbation equations for a 
general model of dark energy and an arbitrary form of coupling to dark matter. We pay special attention to the dark 
energy sound speed and pressure fluctuations. In Sec. Ill we present a simple physically motivated coupling, and we 
derive analytical solutions for the perturbations in the early radiation era, in the case w x = const. These solutions 
reveal the non-adiabatic super-Hubble instability, which is confirmed by numerical solutions using a modified version 
of CAMB Q. In Sec. IV, we extend the analysis to another simple class of background couplings, and we show that 
the same instability persists. We conclude in Sec. V. 

II. DENSITY AND VELOCITY PERTURBATION EQUATIONS 

The general perturbation equations for coupled fluids are given in [10( and various subsequent papers. We follow 
broadly the notation of [ll| and specialize to the case of dark sector coupling. We pay particular attention to the 
covariant form of the coupling and the correct treatment of momentum transfer (which vanishes in the background). 

Scalar perturbations of the flat FRW metric are given in general by 

ds 2 = a 2 1 - (1 + 2<j})dr 2 + 2d l B drdx 1 + (1 - 2ip)6ij + 2d l d J E dx l dx J \. (6) 

The background four- velocity is u^ = a~ 1 8 t Q , and the A-fluid four- velocity is 

u l " A = a- 1 (l-(f>,d i VA), u^ = a^-l-4>,Si[vA + B]j, (7) 

where va is the peculiar velocity potential. The volume expansion rate, which generalizes the Newtonian relation 



= V • v, is [12 1 



6 A = -k 2 (v A + B). (8) 

Energy-momentum tensors 

We choose u A as the energy- frame four- velocity, i.e., there is zero momentum flux relative to u A , so that T Al/ u A = 
-Pau a . Then the A-fluid energy-momentum tensor is 

TH v = ( PA + P A )u^ + PAS> i u + ^ v , (9) 



where pa — PA + $Pa and, Pa = Pa + &Pa ■ The anisotropic stress 7r^ ^ is given by 

1 



^Av = °i 



<9 J <9 



, - -«f , V 2 ) tt a . 



(10) 



'.V, { : 3 

The total (conserved) energy-momentum tensor is T^ v — ^2 T Al/ , so that 

(p + P)u^u v + P^ v + i^ v + g"u„ + q v u» = J2(PA + PA)U* A U* + J] P A + E ^ 

A A A 

Here g M is the total momentum flux relative to the total four- velocity u^ . In general this four- velocity has the form 



Av ■ 



(11) 



u> x = a- l (l-<t>,d i v\. 



(12) 



The choice of v depends on how the total four-velocity is defined. 

It follows from Eqs. (jlip and (fT2f that p = J^PA , P = Yl Pa , ^v — Yl 7r Av > ano - ^ ne total momentum flux is 
q' = a^ 1 J2(pa + Pa)9 1 va — a~ 1 (p + P)d l v. Thus the total energy frame (q 1 — 0) is defined by 



(p + P)v = ^2(pa + Pa)va 
This is the choice of v that we will use from now on. 



(13) 



Energy-momentum balance 
The covariant form of energy-momentum transfer is [10 . 1 131 ] 

v„n»=Q» A , y,Qa = o. 

A 

A general energy-momentum transfer can be split relative to the total four- velocity as [10 . Ill| 

Q A = Q A u^ + F^, Q a = Qa + SQa, V^ = 0, 



(14) 



(15) 



where Qa is the energy density transfer rate and F A is the momentum density transfer rate, relative to w M . Then it 
follows that F A — a _1 (0, 9 1 /a), where f A is a momentum transfer potential, and 



= ad. 



Q A (l + <t>)+8Q A 

f A + Q A (v + B) 



(16) 
(17) 



The perturbed energy transfer includes a metric perturbation term Qa4>> m addition to the perturbation SQa- The 
perturbed momentum transfer is made up of two parts: the momentum transfer potential Qa(v + B) that arises from 
energy transport along the total velocity, and the intrinsic momentum transfer potential j a- In the background, the 
energy-momentum transfer four-vectors have the form 



Q% = a-\Q c ,0) = a-\-Q,0) = -Q£ , 

so that there is no momentum transfer. 

Total energy-momentum conservation implies 

o = E^ = E ,5 ^ = E^- 

For each A-fluid, Eq. (I14p gives the perturbed energy and momentum balance equations (in Fourier space), 
Sp' A + m{8p A + 5P A ) -3(pa + Pa) $ - k 2 (p A + Pa) (v a + E 1 ) = aQ A 4> + aSQ A , 



[{PA + Pa) (va + B)]' + AH (p A + Pa) {v A + B) + {p A + Pa) 4> + SP A - -^tt a = aQ A (v + B) + af A 



2k 2 
3^ 



(18) 

(19) 

(20) 
(21) 



Sound speed and pressure perturbations 

The sound speed c ? a of a fluid or scalar field, labelled by A, is the propagation speed of pressure fluctuations in 



the A rest frame 



eu CsAOi a : 
[Hill HI: 



5Pa 



-sA 



Sp/ 



(22) 



For a scalar field ip, the rest frame is defined by the hypersurfaces ip = const, orthogonal to the rest-frame four-velocity 
m^ ex V^. Thus the kinetic energy density in the rest frame is — ^V^V^i/? = ip' 2 /(2a 2 ), while Sip = in the rest 
frame, so that 5V = 0, where V(ip) is the potential. The density and pressure perturbations are consequently equal 
in the rest frame: 5p v = 5(^a~ 2 ip' 2 + V) = a~ 2 (p'5ip' — 8{\a~ 2 <p' 2 — V) = 8p v . The sound speed is therefore equal 
to the speed of light, independent of the form of V (<£>): 



<„ = i • 



(23) 



We can define the "adiabatic sound speed" for any medium via 



p'a 



= W A 



Pa/ PA 



(24) 



For a barotropic fluid, c 2 = c 2 , and if w = const, then c 2 — w. By contrast, for a scalar field, c 2 ^ c 2 ^ w v . 

The fluid model for dark energy with constant w x is at face value a barotropic adiabatic model. But if we treat the 
dark energy strictly as an adiabatic fluid, then the sound speed c sx would be imaginary (c 2 x — c; 



c ax = w x < 0), leading 
to instabilities in the dark energy. In order to fix this problem, it is necessary to impose c 2 sx > by hand [15| , and it 
is natural to adopt the scalar field value Eq. (|23| . Thus 



1 



w x 



const < . 



(25) 



This is what is done in the CAMB [9j and CMBFAST [JJ] codes. 

In the perturbation equations ([20]) and ([21]) , we need to relate SPa to 6pA via Eq. ([22]) . The A rest frame (the zero 
momentum gauge or comoving orthogonal gauge) is the comoving (wA|rf = 0) orthogonal (-B| r f = 0) frame, so that 



T\ o | r f — 



1 Ai\ ri ■ 



(26) 



We make a gauge transformation, x M — > x^ + 
v a + B = {va+B)\, + Sta, 



(5ta, d l Sx a) , from the rest frame gauge to a general gauge: 



6P A = SP A \ rt -r A 



P'a^ta 



Spa = <W| rf - Pa$ t a ■ 



(27) 



Thus Sta = va + B, and substituting into the pressure and density fluctuations, we obtain 

sp A = c 2 aA 5p A + (c 2 sA ~ c 2 aA ) [6 pa + p'a {va + B)\ = c 2 aA 5p A + 5P nadA , 



(28) 



where 5P n!lt iA is the intrinsic non-adiabatic pressure perturbation in the yl-fluid. 

Our result applies to both coupled and uncoupled fluids, where the difference enters via the term p' A . This recovers 
the expression for the uncoupled case in [lj, Qj|. For the coupled case, the background coupling Qa enters SPa 
explicitly: 



6Pa = 



sA 



Sp/ 



-sA 



L-aA 



STi{l + w A )pA - aQ A 



(29) 



This corrects the expression used in [17| , which omits the Qa term. As a consequence, there are errors in the equations 
in [13, [l8[ for S' x and 9' x . We will discuss the implications of this in Sec. V. 



General equations 



From the above equations we can derive evolution equations for the dimcnsionless density perturbation 8a — Spa/ Pa 
and for the velocity perturbation 9a (which has dimension of k): 



s' A + m(c 2 sA - w A )8 A + (i + w A )6 A + m[m{\ + w A )(c 2 sA - «u) + «4] % 



3(1 + waW + (1 + w A )k 2 (B - E') 



aQA 
PA 



9' A + H(l - 3c 2 sA )9 A - ,° A - k 2 S A 

[1 ■+■ WA) 



3a 2 (l + wa)pa 



8 A + 3H(c 2 A - w a 



k\ A - k 2 



k 2 



k 2 



PA 



■SQj 



aQ A 



(1 + w a )pa 



[e-(i + c 2 A )e A ] 



(1 +w A )p, 



■k 2 f A 



(30) 



(31) 



The curvature perturbations on constant-/^ surfaces and the total curvature perturbation (on constant-p surfaces), 
are given by the gauge-invariant quantities 



U = -v>- ^-—r . t = -v>- K— = > -y (a ■ 
Pa P P 



The total energy conservation equation leads to 



? = -, 



n 



■8R 



ii ad 



(p + P) 

The gauge-invariant relative entropy perturbation for any two fluids is 

Sab=3h( S -^- S -^-)=HCa-Cb). 
\Pb Pa J 



(32) 



(33) 



(34) 



III. A COVARIANT MODEL OF DARK SECTOR COUPLING 



In order to apply the equations of the previous section, we need to choose a model of the dark sector coupling via 
a covariant choice of the transfer four-vector Q% = —Q%- For example, a coupling model motivated by scalar-tensor 
theory has [19| 



-Q^%)T C "„W, 



(35) 



where ip is the scalar field dark energy and j3 is a coupling function. Using this form, the perturbed energy transfer 
8Q C = —8Q X and momentum transfer f c — —f x can be calculated unambiguously. Note that V M <yj is parallel to the 
dark energy four- velocity u£, i.e. Qc — ~Qx ^ u x- 

For more phenomenological models of coupling, especially in the context of fluid dark energy, it is not always clear 
what the covariant form of the transfer should be. For example, consider the background transfer models [201 ] 



n U ( 

Q = —\0-cPc 

a 



a x p x ) 



(36) 



where a c and a x are dimensionless constants. One problem with these models is ambiguity under perturbation: What 
is the covariant form of Q% = —Q^ that reduces to Eq. (|36[) in the background? This has not been made explicit in 
previous work [17], [l8| , as pointed out in [3] . 

A further problem with Eq. (|36[) is the explicit presence of the universal expansion rate Ti. This is designed 
for mathematical simplicity rather than physical motivation: one does not expect the dark sector coupling at each 
event to depend on the global behaviour of the universe, but to depend only on purely local quantities. A generous 
interpretation is that the TL factor is an approximation to the temperature-dependence of the interaction rate. 

Here we propose a covariant model of coupling that avoids these problems: 



-Q£ = r:c< = -rp c u£, 



(37) 



where T is a constant interaction rate, p c is the dark matter density in the inhomogeneous universe, and u£ is the 
dark matter four-velocity. The notable features of this phenomenological coupling model are: 

(1) The interaction rate T is 'local', i.e. it is determined by local interactions and not by the universal expansion rate. 

(2) In the rest frame of the dark matter, there is no momentum transfer. (By contrast, for Eq. (|35[) the momentum 
transfer vanishes in the dark energy rest frame.) 

(3) The case T > corresponds in the dark matter frame to the decay of dark matter into dark energy. This opens 
the possibility of an alternative approach to the coincidence problem: instead of trying to achieve a constant nonzero 
ratio Cl c /Q, x on the basis of primordially existing dark energy, one could try to build models where the dark energy 
accumulates via the decay of dark matter, and dominates in the late universe because the decay rate T is small. 

Background dynamics 

In the background, the coupling ([37)) reduces to Eq. (fl"8[) . with 

Q = Tp c . (38) 

When r > 0, this coincides with a special case of a model in which superheavy dark matter particles decay to a 
quintessence scalar field [21[ . It also has the same form as simple models to describe the decay of dark matter into 
radiation [22], or a curvaton field into radiation [llj. The background dynamics with Eq. (I38|) have been analysed 
in [j] for the case of scalar field dark energy. The growth factor and weak lensing have been investigated for the case 
w x (a) = Wo +w a (l — a) in [23j. 

We impose the condition w x > — 1 so as to avoid a phantom fluid model of dark energy. In order to have a close- 
to-standard matter dominated era for structure formation, and in order to be consistent with the observed angular 
diameter distance to last scattering, it is necessary that | if e ,efif I is small, i.e. 

— <0.1H , (39) 

Pc 

where H = Ho/a, = Ho is today's Hubble rate. 
For this coupling model, Eq. $5§ implies that 

aT aT p c 

W c ,eff=^r7, W Xte $-W x = -— — , (40) 

ort ort p x 

where |T| -C Hq by Eq. ([3"9"|) . Since w x > — 1, it follows that the total effective equation of state satisfies Wtot > — 1> 
so that H is a decreasing function and therefore a\T\/H decreases as we look backward into the past. Thus w c ,<$& <C 1 
for all times up to the present, i.e., the dark matter effectively does not see the coupling for all times from today to 
the past: 

p c = p c0 a- 3 . (41) 

If |r|/-Ho < 3(1 + w x )£l x o/£l c o, then the dark sector coupling is negligible at late times, by Eqs. (|4"0"[) and ([41]) . and 
we have p x = p^o^ 1+ • This is valid to the past until the coupling term |arp c | is equal to the redshift term 
|37Y(1 + w x )p x \ in Eq. ^. At earlier times, the coupling term will dominate for a small enough a, regardless of how 
small |r| is. 

In the radiation era, 

H = t~ 1 , a 2 = HlVLror 2 , (42) 

and the energy balance equations fl} and @ lead to a simple solution for early times in the case w x < —2/3: 

aT^ = (3w x + 2)t- x , w x <-\. (43) 

Px 3 

From now on we assume that w x < —2/3, which is consistent with observations. Equations (|4ip and (|43|) imply that 

2 
p' x = -Hp x =4- W x , e tt=-~ and p x ex a -1 ex t _1 . (44) 



Equation (|41j) shows that p c > for all times. Therefore by Eq. (|43|) , p x becomes negative 1 when Y > 0. This is the 
case that corresponds to the decay of dark matter to dark energy. The rigidity of the assumption that w x is constant 
leads to this problem with the decaying dark matter case. When dark energy is modelled as a scalar field [J], it 
remains positive for all times and there is no such problem in the case Y > 0. 

These analytical approximations are confirmed by numerical integration, as illustrated in Fig. [TJ 

Dark sector perturbations 

For the model described by the covariant energy-momentum transfer four-vector Q£ = — Q£ = —Yp c u^, we need 
to determine SQ C = —SQ X and f c — —f x from the conditions imposed by energy-momentum balance. Using Eqs. ([7|) 
and (|37p. we find the components of Q^ = —Q^: 

Q^ = -Ql=aY Pc [l + <)) + S c , Qiivc + B)]. (45) 

Comparing with Eqs. (HHJ) and (fT7|) . it follows that 

SQ C = -Yp c 5 c = -8Q X , f c = Yp c (v - v c ) = -f x , (46) 



where v is the total energy frame velocity, defined by Eq. (|13p . 

The density and velocity perturbation equations (|30j) and (|31[) for the dark sector, with n c = = n x , w c = = w' x 



and c 2 sx — 1, can then be given, in longitudinal (Newtonian) gauge (B = E = 0): 








5' x + 3H(l-w x )6 x + (l+w x )e x + m 2 (l-w 2 x )%-3(l + w x W = aY^ 

k z p x 


5 C - S x + 3W(1 - 


■w x )-^ + cf> 


. (47) 


k 2 nY n 

e' x -2H9 x - S x -k 2 cl> = £l(e c ~29 x ) 7 
(l+w x ) (l + w x )p x 


(48) 


and 




5' C + C -W = -aY4>, 


(49) 


e' r +ne c -k 2 d> = o. 






(50) 



Note that the dark matter velocity perturbation equation (|50p is the same as in the uncoupled case. (In particular, 
this means that in synchronous gauge, we can consistently set 9 C — 0, as is done in the standard, uncoupled case.) 
This is due to the fact that there is no momentum transfer in the dark matter frame. 

We are now in a position to see qualitatively why there is a large-scale instability in the dark sector perturbations 
during the early radiation era. The coupling term Q x in 5P X , Eq. (|29[) . leads to a driving term 

„ aY p CQ -2(3w x + 2)_ mx 



(l + w x )p x l + w x 

on the right hand side of Eq. (|48p. Here the multiplier oiTL9 x is a positive number (since w x < —2/3) - and it becomes 
very large if w x is close to — 1. This causes rapid growth of 9 X . Qualitatively, this is the source of the instability: 
in the presence of energy-momentum transfer in the perturbed dark fluids, momentum balance requires a run-away 
growth of the dark energy velocity. The precise form of the instability is computed analytically below. 

Radiation era 

Tight coupling between photons and baryons means that (a) the only nonzero momentum transfer is in the dark 
sector, and (b) the only nonzero anisotropic stress is that of the neutrinos (which have decoupled). The perturbed 



1 Note that we can avoid p x < for T > if w x > —2/3 or w x < — 1, but we exclude these cases. 



Einstein equations reduce to 

3t _ V + k 2 ij) + 3T- 2 (f> = -4nGa 2 5p, (52) 

/fc 2 (^' + r-V) = 4irGa 2 (p + P)6, (53) 

k 2 

■ijj" + 2r _ V -t^iP + t- 1 ^ + —(i/j-<f>) = 4irGa 2 6P, (54) 

tp-cf) = 87rG7r„. (55) 

The perturbed balance equations in the dark sector are given by Eqs. (|47|) - (|50p . with background coefficients deter- 
mined by Eq. (|42[) . For the photon-baryon sector, 

5' = --0+4a/)', 8' b = -6 b + ty, (56) 



3 "> 
f 

4 



■= -fc 2 £ 7 + fc 2 </>, 9' b = -H6 b + c 2 b k 2 5 b + k 2 4>, (57) 



and for neutrinos [li 



<^ = --0„ + 4<//, e' v = - A k 2 5 v + k 2 <i>-k 2 o- v , <7t = 4 fl «" ( 58 ) 

3 4 15 

where o v :— 2k 2 "K v /[3a 2 (p v + Pi/)], and we have neglected the neutrino octopole, i.e., we work to leading order in kr. 

Adiabatic initial conditions 

Now we look for a solution in the radiation era, in the super-Hubble scale limit, kr -C 1. We find that we can set 
adiabatic initial conditions to lowest order in kr: 

4> = A 4> = const , ij> = ( 1 + ?rJ) <t> , R»~ PV , (59) 

V 5 ) p u + p~f 



8^ = 5 U = -5 b = -s c = -20 , e 7 = e u = o b = e c = -(kT)k(t>, G V = —{krf^ 1 (eo) 

6 X = i(5 7 , 6 X = 7 . (61) 

The expression for 5^ follows from £ x — C 7 = [see Eq. ([34]) ] , using the early-time attractor solution Eq. (|44| . 

However, at higher order in fcr, this solution is not adiabatic. In the standard case with uncoupled dark energy, 



this is not an issue - since the deviation from adiabaticity is decaying and suppressed 24|. However, for the coupled 
dark energy model considered here, the situation is dramatically different - because there is a strongly growing non- 
adiabatic mode on super-Hubble scales. Even with adiabatic conditions in the limit fcr — > 0, the terms of higher order 
in kr contain the non-adiabatic mode, and this mode will dominate since it is strongly growing (see below). Note 
that this mode is regular, and it is stimulated by the dark sector coupling, as explained via Eq. (|51|) . 

The detailed analysis of all perturbative modes is given elsewhere [25J . There we also show numerically how the total 
curvature perturbation starts off constant with the initial conditions Eqs. (f59"|) -([6"i" |) . but begins to grow dramatically 
after a short time (well before equality). 

The dominant non-adiabatic mode 

In order to find the non-adiabatic mode, we assume a leading-order power-law form for the perturbations: 

$ = A+Qst)"* , ^ = A 4> {krY^ 6 A = B A (kT) n \ 6 A = C A {k T y A , a v = D v {kr) n ° , (62) 



where n^p and n$ are not zero (n^, = = n<f, is the adiabatic case). The perturbed Einstein and balance equations ([47)) 
(|58p may then be solved, to leading order in fcr, in terms of i\>: 

<j) = Jip, (63) 

S y =5 u = 4a/), (64) 



5 1 = 3ip, (65) 

(7 " : " ' l '-(kr)ktl>, (66) 



>v — Vb = 



(fty + 1) 



' (*r) % , (67) 



15(n^ +2) fc 

{jH + !) J g 

(n^ + 2) (J + l) 7 ' 

2fc 3 (3w c - 
(n^ + 2) 



(68) 



2k 3 (3w x + 2) [n| + (J+1K - (J + 2)] , 



«" = rfe K-V (fcT) "*■ (69) 



(kr)kd x . (70) 

3(1 -Wj) 

The gauge-invariant curvature perturbation, Eq. (132[) . is given in terms of ip as 

C = -~(n^ + J + 2)V>. (71) 

The solution is thus fully determined up to an arbitrary normalization of the amplitude parameter A^. The stress 
anisotropy parameter J is given by 

J .= i* = 1 1^ (72) 

A^ 5(n v , + 2)(n 1/J + l) + 8i?„ 



The power-law index n^ is determined in terms of w x as 



~(l + 2 W:c )±73^^2 

"0 = n± = — . (73) 

1 + w x 

The fastest growing mode is the n + -mode. Equation (f69|) shows that the modes are regular (i.e., well-behaved as 
kr — > 0) provided that Ken± > 3. This leads to the conditions, 

4 .9 4 

n + regular if — 1 < w x < — , n_ regular if < w x < — (74) 

5 11 5 

The only explicit dependence on the coupling rate T in the early-time solutions is in the 5 X solution, Eq. (|69[) . The 
uncoupled limit Y — leads to a singularity in that equation. This reflects the fact that the solutions are not valid in 
the limit Y = 0. We cannot recover the r = limit since we have used in an essential way that r/0, see Eqs. (|43|) 
and (|44[) . The most important difference is that n^ = in the uncoupled case, so that tp and (f> are constant, and in 
addition £ = const. By contrast, when Y ^ 0, we see that n + is typically very large: 

2 , , 

w x 1 => n + - > 1 . (75) 

1 + w x 

This large exponent signals the blow-up of ip, and therefore of all perturbations, including the gauge-invariant curvature 
perturbation, Eq. (|71|) . on super-Hubble scales in the early radiation era. The instability is stronger the closer w x is 
to -1. 

The instability occurs no matter how weak the coupling is. A smaller value oflrl simply moves the blow-up to 
earlier times. This is in contrast with the strong-coupling instabilities discussed in [7J,l8|. Furthermore, the instability 
is non-adiabatic, in accordance with Eq. ([3"3"|). since the curvature perturbation blows up. Again, this is in contrast to 
the adiabatic instabilities of 0, Q . The origin of this large-scale non-adiabatic instability is not simply the fact that 
the dark energy fluid is non-adiabatic, i.e., (? ax ^ c 2 sx . In the uncoupled case, the same non-adiabatic fluid behaviour 
is also present, but there is no instability. The coupling plays an essential role in driving the large-scale non-adiabatic 
instability. 
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Coupling Q = Tp c with I7H = -5e-7 



Coupling Q = Tp with w = -0.87 
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FIG. 1: The background evolution of dark energy density p x (top panels) and the evolution of the gauge-invariant curvature 
perturbation £ for a super-Hubble scale k — 7 x 10~ 5 Mpc -1 (bottom panels), as functions of scale factor a, for the model 
with coupling given by Eq. (|37|l . In panels on the left, we vary w x with F fixed, while the right-hand panels have fixed w x and 
varying V. For each case, the vertical lines (in matching style) indicate the moment when |37i(l + Wx)px\ and jarp c | are equal 



in Eq. ([2|. To the right of these lines the background evolves as in the uncoupled case, i.e, p x ex a 



-3(l+u> x ) 



To the left, the 



coupling modifies the background evolution to p x oc a -1 , Eq. (|44[) . (Note that in the left panels for w x = —0.80 this happens 
very far in the past, not shown in the figure.) All the curves show the full numerical solution obtained with our modified version 
of CAMB, with the initial amplitude of if) set to 10 -25 . The analytical solution for the blow-up of £, Eq. (|71|l . is practically 
indistinguishable from the numerical solution at early times. 



The analytical demonstration of the large-scale instability at early times is confirmed by numerical solutions, using a 
modified version of CAMB. Examples are shown in Fig.[TJ As long as the coupling modifies the background evolution, 
the curvature perturbation is extremely rapidly growing: Q ex a n+ where n+ is given by Eq. (|73p during radiation 
domination. When the background starts to behave as uncoupled, the blow-up of C ends and it begins to oscillate 
about zero with large (and mildly) increasing amplitude. 

Equation (|73[) shows that for large enough w x , there are oscillations super-imposed on the super-Hubble blow-up 
mode: —y/2/3 < w x < —4/5, where the upper limit comes from Eq. ([74]) . See Fig. Q]for an example (left panel, with 
w x = -0.8). 

Figure Q] also shows that if |T| is large enough to modify the background after radiation-matter equality, then |£| 
follows a matter-dominated attractor solution: 



C oc a n+ where h + ~ — n+ for w x <- 
The background solution in the matter era that corresponds to Eq. (|43l) is 



1. 



(76) 



aT— = 3(2^ + ly- 1 

Px 



w x <-- 



(77) 



The scale-dependence of the instability is illustrated in Fig. [21 The important point is that the instability cannot 
be removed by simply re-scaling Aj, : even if we can match the large-scale CMB power with a small enough A^ , the 
full CMB and matter power spectra will exhibit a strong scale-dependence in violation of observations. 
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Coupling Q = rp with r/H = -5e-7 and w = -0.87 



— k = 7e-5 Mpc 
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k = 5Mpo" 1 



Coupling Q - Tp c with r/H = -0.1 and 



-0.87 





FIG. 2: The evolution of the gauge-invariant curvature perturbation £ for three different scales as a function of scale factor 
a, for the model with coupling given by Eq. (|37p . In the panel on the left, the coupling F| is very small while in the right- 
hand panel, F| is larger. Vertical lines indicate the moment when each mode enters the horizon (fcr ~ 1). The largest scale 
(k = 7 x 10 -5 Mpc -1 ) stays super-Hubble all the way up to today. The intermediate scale (k = 1.5 x 10~ Mpc -1 ) enters the 
horizon during matter domination, and the smallest scale (k = 5 Mpc -1 ) enters deep in the radiation era. 



Note also that the analytical derivation of the instability does not depend on the sign of T, i.e., there is a blow-up 
for positive and negative I\ However, as shown by Eq. (|42[) . when T > 0, there is always a time r* when p x goes 
through zero and is negative for t < t*. The perturbation equations are singular at r*. 



IV. EXTENSION TO OTHER COUPLING MODELS 



The instability in the model with coupling Q^ = — Q£ = —Tp c u£ is not peculiar to the particular form of the 
coupling. The other models discussed in Sec. Ill also suffer from this instability. We consider two background 
couplings, which are special cases of Eq. ([36]): aQ = aHp c and aQ = /3H.(p c + p x ). 



Model with background coupling aQ — aHp c 



For this coupling the background balance equations may be solved exactly: 

,-(3+a) 

a 



Pc = PcQO, 

Px = Px oa- 3{1+w ^ + 



a — 3w a 



p c0 a (a 



-3w x 



(78) 
(79) 



where we assume that a > 3w x (otherwise the coupling strength |a| would be too large). The dark matter density 
is always positive. For the dark energy density to cross zero and become negative, there must be a solution a n < 1 
to a^""~ a = 1 + Q x n{a — 3w x )/(Q c oa). The left hand side is always > 1, since 3w x — a < 0, whereas the right hand 
side is > 1 only for a > 0. Thus for a > 0, i.e. for the case of dark matter decaying into dark energy, the dark energy 
density always becomes negative in the past. In the early universe, a <C 1, the exact solution implies 



Px 

Pc 



3w x — a 



a > 3w x . 



(80) 



To analyze perturbations in this model, we need a covariant form of energy-momentum transfer four-vector that 
reduces to Q — aTip c /a in the background. We propose to use the same form as Eq. (j37|) : 



'dH.p c u p c 



(81) 



Here we are following the implicit assumption made in other work that aH is an approximation to an interaction rate 
that varies with time but not in space, so that there is no perturbation of H in 6Q^. 
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With this covariant form of energy-momentum transfer, the momentum transfer is 

af c = aHp c (v - v c ) = -af x , (82) 

and the dark sector density and velocity perturbations are given by Eqs. (|47 p - (|50j) . with T replaced by atija. 

The early radiation solution to leading order in kr is qualitatively similar to the solution for the T model in the 
previous section, with differences arising because the interaction rate aH varies with time, as opposed to the constant 
rate T. The solutions for the dark energy perturbations become: 

2v^(3wg-a) [n% + (J+l)n+-(J + 2)] k 2 

®x = q n n - — v — -, — — n 7T^' (83) 

3aU c0 {l + w x ) (n + + l) H 

^ = - 3( r + i )( 2 1 "? ) ^" i T- (84) 

(n+ + 2 — a) k 

The key indicator of instability, i.e. the power-law index for the fastest growing mode, n+, takes a more complicated 
form: 



3a + (a - 6)w x + J (a 2 + 16a + 40) w;^ - 2(a 2 + 6a + 8)w x + (a 2 - 4a - 20) . , 

n + = 7T, s ■ (85) 

+ 2(1 + w x ) v ' 

For example, if w x = —0.87 and (3 = —0.003, then n + = 38.95. Unlike the T-model expression Eq. ([75)) . here n + 
depends explicitly on the interaction rate parameter a. Note that, as for the T-model, the limit a = is not admitted 
in Eqs. (|83|) - ([85)) . since the derivation uses a ^ in an essential way. For small \a\ and w x close to — 1, 

f ' » 1 . (86) 



1 + w. 



This is triple the corresponding index for the V model. The analytical form for the early-time instability is confirmed 
by numerical integration. 

Model with background coupling aQ = f57i{p c + p x ) 

The background coupling Q is proportional to the total dark sector density p c + pxi which obeys an energy conser- 
vation equation. The balance equations lead to an exact solution [20j for p x /p C7 with 3 cases according to the sign of 
— (3 + 3w x /4. The non-negative cases are not relevant since they violate |/3| <§C 1. For the remaining case, 

p x (B + 2/3 - 3w x ) ( 1 - ba B+ W- 3w * \ , 3w x - 2/3 3w x 



p c 2/3 \l + ba B + 2 f 3 - Sw * J 2/3 4 

where 



B := ^3w x {3w x - 4/3) - 2/3 + 3w x 
and b := [B - 2[3Q, x0 /n c0 ]/[B - 6w x + 2/3(2 + n x0 /£L c0 )]. It follows that 

Px 1 



p c 2/3 X { B - 2^3^(3^-4/3) a -* oo (89) 

If /3 > 0, then p x /p c is negative in the early universe, a —> 0, and also in the future, a — > oo. Therefore the /3 > case 
of this model is unphysical. For /3 < 0, if we fix f^o/^co a t a value greater than the late attractor in Eq. (|89")1 , then 
in the past p x j p c becomes negative. For a physical model, we thus require 3w x /4 < /3 < and Q x o/fl c o less than the 
late-time attractor. 

Perturbation of this coupling model is more complicated because it is determined by the total density, and therefore 
there is more ambiguity in the appropriate choice of four- velocity in the definition of Q^ — —Q x - In previous 
work [13j[3' this issue was not explicitly discussed, and no form for Q% = —Q^ was given (this was also pointed out 
in [7|). It appears that the dark sector perturbation equations in [17| do not conform to momentum balance. They 
also neglect the coupling term in the expression for SP X , i.e., the Q x term in Eq. (|29|) . as we pointed out in Sec. II. 
It turns out that this second error is decisive for the instability, whereas the error in f c = —f x only leads to small 
corrections. 
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Using the correct form, Eq. (j29f . for 5P X , the evolution equation for dark energy density perturbations (which is 
independent of f x ) becomes 



S' x + 3W(1 - w x )5 x + (1 + w x )0 x + m 2 {l -wl)-^- 3(1 + w x W 



■■0H 



<fi + m(l-w x )-%} + ^(5 c -6 x ) 

Px J L K J Pa; 



(90) 



In [13, [3, Eq. (7), the right hand side has — S c — p x 6 x /p c instead of S c — S x . We find an instability (see below), 



whereas there is no instability in the results of [lTf ■ Their omission of the coupling term in 8P X has inadvertently 
removed the instability. 

The dark matter density perturbations obey 



5' +e c - 3v/ = -m 



1 



fh_ 

Pe 



■ + ^(S X -S C ) 

Pc 



(91) 



The equations for 5 X and 5' c are independent of the momentum transfer f c — —f x . In order to compute the velocity 
perturbations, we need the momentum transfer. If we follow the previous models and choose the energy-momentum 
transfer four- vector to be aligned with the dark matter four-velocity (so that there is no momentum transfer in the 
dark matter frame), then 



aQ£ = -aQZ = -PH{pc + p x )u» 
It follows that 

af c = (3H(p c + p x )(v - v c ) = -af. 
and the velocity perturbation equations become 

k 2 „ (iU 



2H6 X 



-S T . — k d> = 



(l + O x 
e' r + ne c -k 2 (h = 0. 



l + ZL)[B e - 20 x ] , 
(1 + w x ) \ p x • 



(92) 
(93) 

(94) 
(95) 



The 0' c equation agrees with |17l.[l8|. but their 9' x equation has +9 X in place of our 9 C — 28 x on the right-hand side. 
Any other choice for the four- velocity along Q£ will lead to a nonzero right-hand side in Eq. (f95|) . It appears that no 
consistent choice of Q 1 ^ can recover the equations of [13, [3 ■ 

We find that the problem of not accounting correctly for momentum transfer has a minor effect on the instability. 
The key driver for the instability is the coupling term in SP X , leading to the correct forms Eq. (|9"0"1) and ([M| for the 
8 X and 8' x equations. 

Using Eqs. (|9"0|) - (f9"5)) . we can find the early radiation solution. The solution is simplified if \/3\ is small enough, and 
we find that: 



4/VTU \n% + (J + l)n+ - (J + 2)1 k 2 , 

-no 



wn c0 (i + w x ) 

(1 + w x )(n+ - 2) + 



(n+ + 1) 
2/3(B + 2/3) 
B 



M'T 



(96) 

(97) 



The power-law index is given by 
1 



" + 



2{l + w x ) 



M(l+w x ) -2N 



+ \ M (1 + w x ) + 2N - 4 (1 + w x ) (M + 2) (-2 - 2w x + 2N) - 3 (1 - w x ) (N - 3 - 3^) 



(98) 



where M = -3w x + 2/3 2 /B and N = (3 + 2f3 2 /B. For example, if w x 
in the previous model. For \j3\ <SC 1 and w x ~ — 1, we have 



6 

1 + W a 



>1 



= -0.87 and = -0.003, then n+ = 38.95, as 



(99) 
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Coupling Q = p H (p +p x ) with p = -3e-3 




FIG. 3: The evolution of the background p x /p c (top panel), and gauge-invariant curvature perturbation £ for a super-Hubble 
scale k = 7 x 10 -5 Mpc -1 (bottom panel) for the model with coupling given by Eq. (|92p . The figure shows a full numerical 
solution for w x = —0.87 and /3 = —0.003, starting with the initial value ( = 10 -80 and ending up with oscillations of amplitude 

ici > io +3o °. 

Hence our analytical solution again shows an instability, and this is confirmed by numerical solution, as illustrated in 

Fig.H 

We have also checked numerically that the instability persists with negligible changes if we make alternative choices 
of Qg = —Q%. The two obvious choices are to align Qg — — Q£ along the dark energy four- velocity or along the 
"centre-of-mass" four- velocity in the dark sector: 



-aQ" x = -pn{ Pc + PzK . a,Qc = -aQ'i = -pH(j> e + P x )< x , 



(100) 



where 



(Pc + fe)< = PcK + PxK ■ 



(101) 



CONCLUSIONS 



We have given a detailed general analysis of the relativistic perturbations for a cosmology with coupled dark 
energy and dark matter fluids, paying particular attention to the non-adiabatic features in the dark energy sound 
speed and to the correct three-momentum transfer in the dark sector. We specialized to the case of a constant dark 
energy equation of state w x , and with energy-momentum transfer four-vector of the form Q^ = — Q^. = —Tp c u^, 
where T is the constant interaction rate. We were able to find the fastest growing early-radiation regular solution of 
the perturbation equations to leading order in fcr, and this solution shows a strong blow-up of the gauge-invariant 
curvature perturbation on super-Hubble scales. Even if adiabatic initial conditions are chosen, the dominant mode 
quickly overwhelms the adiabatic mode. 

The instability occurs no matter how weak the coupling is. Decreasing |r| only shifts the blow-up to earlier times. 
This is in contrast with the strong-coupling instabilities discussed in [7|, |8j. The instability is also non-adiabatic, 
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since growth of |£| is driven by the total non-adiabatic pressure. This is also in contrast to the adiabatic instabilities 
of@,|. 

The origin of this large-scale non-adiabatic instability is not simply the fact that the dark energy fluid is non- 
adiabatic, i.e., c„ jt c 2 sx . In the uncoupled case, the same non-adiabatic fluid behaviour is also present, but there is 
no instability [15| . The coupling plays an essential role in the large-scale non-adiabatic instability. It appears that 
the key driver for the instability is the coupling term that enters the non-adiabatic dark energy pressure perturbation 
5P X . This leads to a runaway growth of the dark energy velocity, in order to maintain momentum balance in the 
presence of energy-momentum transfer between the perturbed dark fluids. 

We also showed that the instability is not specific to the coupling model that we introduced. Similar coupling 
models show the same qualitative behaviour. We proposed and analyzed covariant forms of Q£ = — Q% that reduce in 
the background to two previously studied coupling models, aQ — aTLp c and aQ = f3TL{p c + p x ), and showed that in 
both cases the non-adiabatic large-scale instability is present. The perturbations in the second case were previously 
considered in [13, [l8[ , but they omitted the coupling term in SP X and this inadvertently removes the instability. We 
confirmed numerically that different choices of Q£ = — Q%, i.e., aligning them along different four- velocities, has a 
negligible effect on the blow-up of the perturbations. 

Uncoupled models with constant w x and c 2 . x = 1 are perfectly well behaved. But it appears that these models are 
unstable to the inclusion of coupling, at least for simple forms of coupling. What are the ways to avoid this instability? 
We can relax the assumption w' x = 0, using a quintessence model instead of a fluid model for dark energy. Or more 
simply, we can use the parametrization w x = Wq + (1 — a)w a . In this parametrization, w x is effectively constant, 
w x — wq + w a , in the radiation and early matter eras. Therefore, our analysis applies unmodified also to this case. 
This may place a tight theoretical lower bound on w a , since Wq is the late-time value of w x . Whether we can produce 
good fits to CMB and matter power data is under investigation [25[ . 

The results of this paper, together with previous results on adiabatic strong-coupling instabilities, strongly con- 
strain the model space for coupled dark energy. Further constraints on certain models arise from the background 
dynamics [3] and from "fifth- force" - typ e limits in cases where couplings extend to non-dark particles (e.g. in the case 
of supergravity-based quintessence [26|). 
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